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ABSTRACT 



Total heat transfer coefficients were determined for a heated, 
copper cylinder placed horizontally in a fluidized bed. The fluidized 
particles were glass beads ranging in size from 215 pm to 3. 40 mm. 
Fluidization was accomplished using air superficial velocities from 
0.045 kg/m^-s to 2.355 kg/m^-S. Temperatures were measured using 
thermocouples and the resulting h values tabulated and plotted as a 
function of superficial velocity and relative velocity. Nusselt 
numbers based on particle diameter and tube wall diameter were also 
determined. The results show consistent patterns with other workers 
and are substantiated by empirical correlations. 
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INTRODUCTION 



Fluidized Bed Combustion is currently being studied intensely due 
to its very promising potential for producing clean and relatively 
inexpensive power from our nation’s vast coal resources. Beds 
containing limestone or dolomite will absorb sulfur dioxide, thereby 
providing a sensible way to burn high-sulfur coals to generate power 
while maintaining high air quality standards without the use of 
expensive wet scrubbing equipment. Furthermore, the long-term goal of 
a pressurized fluidized-bed coal combustor which will generate steam 
and compressed hot gases to drive both steam and gas turbines in a 
combined power cycle holds the great promise of increasing the plant 
thermal efficiency by fully several percentage points over 
conventional power cycles. 

Heat transfer between the combusting bed and the working fluid 
ordinarily takes place by means of submerged horizontal boiler tubes. 
Knowledge of the magnitudes of the heat transfer coefficients and 
temperature distributions along the tube surfaces is essential for an 
optimum design. The rate of heat transfer to a submerged tube depends 
upon a number of complex factors, including the properties of the bed 
material and the fluidizing gas, bed and tube geometries, and the 
fluidized state [1]. Many investigators have reported measurements of 
overall heat transfer between fluidized beds and horizontal tubes for 
a wide range of operating parameters [2-8]. 

Using a model of bed-to-surface heat transfer in which the total 
heat transfer consists of the additive contributions of gas convection 
and particle circulation, the heat transfer behavior may be determined. 
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Since particle circulation is difficult to measure, its contribution 
can be inferred from knowledge of the gas convective part and the 
total heat transfer coefficient. The gas convective component was 
measured by Matt Wade [9] in a separate set of experiments using a 
mass transfer analogy technique. The technique employs the use of a 
sublimating surface to indicate the convective component. The total 
heat transfer coefficients are determined in this paper and will be 
used for comparison to the gas convective results of [9]. 

Briefly, the investigation involves placement of a hollow, copper 
cylinder horizontally in a fluidized medium of glass beads. The 
cylinder represents a typical boiler tube and the glass beads are 
analogous to the boiler bed. Bed sizes were selected to model typical 
coal particle sizes. The cylinder has a heating cartridge placed in 
the hollowed portion which applies a uniform heat flux. Thermocouples 
are used to measure the bed temperature and resulting cylinder 
temperature. Knowing the heat applied and dimensions of the apparatus 
allow the total heat transfer coefficient to be calculated. 

The results of these heat transfer experiments will be compared 
to the work of several previous investigators [2,4,7]. In addition, 
the present results will be utilized in determining the ratio of the 
gas convective component of heat transfer to the total heat transfer. 
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THE EXPERIMENTAL APPARATUS 



General Arrangement 



The experimental apparatus consists of a compressed air system, 
the containment bed, fluidized particles, the test cylinder and the 
instrumentation package. The air system produces, distributes, 
regulates, monitors and contains the working fluid. A general 
arrangement schematic is presented in Figure 1. Briefly, the figure 
depicts the oil-free compressor as the air source. This compressor 
delivers nearly 350 SCFM at 100 psia. This maximum flow rate was 
sufficient to provide fluidization over the entire range of particle 
sizes used in the experiments. 

The distribution system is composed of 2 in. PVC pipe and 
fittings arranged in a convenient air-tight geometry. Pressure 
regulators, pressure gauges and overall flow meters were used to 
monitor and control the flow. The air then flows through one of two 
rotameters for flow rate measuring. The rotameter selection depends 
on the particle size and associated fluidizing velocity under study. 

Flow control is achieved through the use of gate valves forward 
of the rotameter. The gate valve was selected instead of a globe 
valve to provide greater sensitivity in controlling the flow rate. 
After leaving the rotameter, the flow enters the diffuser where it is 
expanded to the dimensions of the bed itself. The pressure drop 
across the fluidized bed is observed by examining a mercury manometer 
tapped into the entrance of the diffuser. An air distributor system 
spreads the flow evenly to provide uniform fluidization. After 
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passing through the beads, the air is exhausted to the outside by 
means of galvanized duct work. 

Fluidized Bed 



Fluidization is accomplished in an 11 -inch inside diameter 
plexiglass containment cylinder mounted vertically on supports. 

Figure 2 shows the bed in detail. The bed has a 6.5 inch x 9 inch 
removable access window with a matching leak-proof cover. The window 
allows the bed to be filled with particles as well as access to the 
test section. 

The bead fill level is also indicated. This level is a constant 
regardless of the bead size under study. Two drain spouts allow the 
particles to be removed to independent levels. The upper level allows 
clearing only the test section while the lower spout is used to clear 
the entire containment of particles. The spouts are sealed with duct 
tape attached flush to the bed wall. This taping method minimizes 
flow disruption. The flow is left undisturbed by using smooth, 
polished inner walls and by having any protrusions away from the test 
section. 

The base of the bed has a conical, 2 inch to 11 inch sheet metal 
diffuser acting to expand the compressed air. The length provides 
gradual, turbulent-free expansion. A distributor plate made up of a 
packed bed of 3.^0 mm diameter glass beads sandwiched between two fine 
mesh brass screens forms the air distributor between the bed and the 
diffuser. The screens are 3 inches apart and secured by tack welds. 
The mesh of the screen is small enough to prevent the smallest 
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particle size, P-010 from falling through it into the air distribution 
system below. 

The interface between the diffuser/upper screen and the bed is 
made by a gasket-filled flange connection. The gasket is a composite 
of 1/16 inch rubber material layers and tube type sealants. The 
flanged gasket seal maintains the air-tight integrity of the bed. The 
actual bed is shown in Photo 1 . The advantage of using a plexiglas 
containment is apparent from the photograph. Namely, the extent of 
fluidization and general fluid behavior is clearly visible. Photo 2 
also shows the bed when fluidized. 

Fluid Particles and Rotameters 



The fluidized particles are industrial quality glas spheres of 
various sizes. The glass is made inert from its lead-free, soda-lime 
silica composition. The spheres have uniform roundness and each size 
falls within a narrow band of dimensional variation. Since larger 
particles require successively greater volume flow rates to fluidize 
them, rotameters of different flow ranges must be used to maintain 
accuracy. Three meters were used having ranges of 1.35 SCFM to 27 
SCFM, 12.82 SCFM to 1 60 SCFM and 100.5 SCFM to 1005 SCFM. They will 
be identified as small, large and blue, respectively. The first two 
occupy interchangeable positions and are shown in Photo 3 and Photo 4. 
The blue meter is identified by its color and is the monitoring meter 
in the distribution system. Photo 5 shows the blue meter. This meter 
is used because its range exceeds that of the large meter and because 
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meters of the same type as the large and small ones are not available 
in the appropriate size. 

Table 1 summarized the bead identifier, average particle size and 
the meter used. 

TABLE 1 : Summary of Particle Name, Bead Sizes and Meter Used to 

Measure Flow 



BEAD 

IDENTIFIER 


AVERAGE 
SIZE (m) 


METER 


P-010 


2.15x10“ 4 


Small 


P-0230 


5.125x10 _J4 


Large 


P-0470 


9.45x10“^ 


Large 

Blue 


A-205 


2.03x10"3 


Blue 


A-340 


3.40x10-3 


Blue 



The readings from the three meters were all converted to volume flow 
rates. In this way, the scale differences of the meters were 
eliminated and comparisons were made on a uniform basis. To better 
visualize the particle sizes, P-010 is about the size of table salt, 
P-0470 resembles large grain sand and A-340 is approximately the size 
of "BB's" for a pellet gun. 

The Test Cylinder 



A one-inch diameter copper cylinder is the focus of the 



experiment. Its seven inch length and 3/16 inch wall thickness were 
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designed to meet the heat transfer requirements of the experiment. 
Figure 3 and Figure 1* show the cylinder's design in detail. Copper 
was selected because of its high thermal conductivity, thus providing 
a uniform outer wall temperature when heated. To measure the 
temperature achieved and to verify the absence of circumferential 
differences, four chromel-alumel thermocouples were attached at 90° 
positions around the circumference. The thermocouple beads are 
attached below the wall surface and the wires run through machined 
grooves. The grooves are filled with thermal epoxy and sanded flush 
to the surface. In this manner, not only are the temperature 
characteristics unaltered but also the flow around the cylinder is 
left undisturbed. 

An electrical heater cartridge is placed inside the cylinder's 
hollow core. The cartridge acts as a radially symmetric heat source 
to generate the desired temperature gradient between the test cylinder 
and the bed. The cartridge matches the cylinder in both length and 
diameter. To ensure effective conductive contact between the 
cartridge and the cylinder, thermally conductive paste was applied to 
the cartridge before inserting it into the cylinder. The combined 
effects of the planned tight fit and the paste minimized any possible 
contact resistance. 

The cylinder is supported on two nylon end caps fastened into 
plexiglass clamps on the bed wall. Figure 5 shows the end caps or 
spacers. Paramount in their design is the matching with the copper 
cylinder of outer diameters, the fastening with recessed set screws, 
their length and the material selection. In this manner, the fluid 
flow is undisturbed and the test cylinder is positioned in the bed's 
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center. Heat transfer end effects may be ignored because of the 
extremely low thermal conductivity of the nylon material. As a 
further precaution, the thermocouple wires and cartridge power card 
run through a hollow spacer and exit the bed directly. 

The assembled unit is shown in Photos 6 and 7. Photos 6 and 7 
show the cylinder removed and partially disassembled for clarity. 

Instrumentation Package 



The instrument package is relatively simple. It consists of a 
variac power supply, a radial, multi-terminal switchbox, a digital 
voltmeter, chromel-alumel thermocouples and an insulated icebath. The 
circuit diagram depicting the component relationships is shown in 
Figure 6. The actual equipment used is presented on Photo 8. 

The theory of operation combines the instruments to extract the 
temperature of the cylinder, the bed temperature and the power 
supplied. Other meaningful values may then be refined from the data. 
Briefly, the variac provides an adjustable current source to the 
heater cartridge. This source is indirectly measured by using the 
voltmeter. The voltage coupled with calibration information about the 
cartridge's resistance produces the power supplied to warm the 
cylinder. The digital voltmeter is a convenient and accurate 
measuring device for both AC and DC voltages present in the experiment. 
The DC potential differences arise from the use of thermocouples. The 
thermocouples use .010 m chromium wire and .010 m aluminum wire with 
teflon insulation. The thermocouples were manufactured in an inert 
gas environment by electric arc fusing. Manufacture using this 
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technique assures measured values corresponding closely to 
standardized Type K temperature conversion tables. 

The radial switchbox is at the heart of the instrument package. 

It provides a common junction for all the thermocouples as well as 
providing easy selection of voltmeter inputs. Insulation packed in 
the switchbox minimizes false DC readings due to ambient temperature 
fluctuations. An insulated icebath is used as a reference for 
thermocouple voltage to temperature conversion. Careful use of 
distilled water and ice provided uniform reference from one data run 
to another. In addition to the cylinder and icebath thermocouples, 
three thermocouples indicate the inside bed temperatures. The 
thermocouples penetrate the bed wall at positions of 1/3 the bed 
circumference and 1/3 the elevation to the test section. This helical 
placement pattern averages out local temperature gradients to yield a 
good overall impression of bed temperature. 

Data collection using the instruments discussed requires manual 
recording of voltmeter readings. Automated "Dataloggers" on the other 
hand, electronically scan, average and convert the readings to 
temperature. Attempts were made using a Monitor Labs 9302 datalogger 
but the results proved unsuccessful. The reason for the failure is 
attributed to excessive static discharge. The static is generated 
from the plexiglass bed's natural interaction with the glass beads. 

The charge then travels down the thermocouple wires into the 
datalogger causing problems ranging in severity from erroneous 
readings to complete microprocessor failure. Mitigating efforts 
included: grounding the bed, grounding the datalogger, installing 

terminal connection slunts, electrically insulating the test cylinder 
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with 4 coats of polyurethane and even applying commercial static 
minimizer sprays. Failure of all attempts to reduce the static to 
within the datalogger's acceptable limits forced the abandonment of 
automated data collection. The digital voltmeter proved to be less 
sensitive to the static and was selected as the logical alternative. 
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FI GURS 1 : Schematic of Overall Flow and Fluidized 

3ed System 




FIGURE 2 s Detailed Schematic of the Fluidized Bed 
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PHOTO 1 



luidized Bed Apparatus 
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PHOTO 2 : Close-up of Fluidized Eed Containment 

Vessel (Test Section) 
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PHOTO 3 : Small Rotameter 




PHOTO 4 : Large Rotameter 
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PHOTO 5: Blue Rotameter 




PHOTO 6 : Test Cylinder Assembled 
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FHOTO 7 : Test Cylinder Separated To Show 

Heater Cartridge 




PHOTO 8 : Instrument Package 
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THE EXPERIMENTS 



Procedure and Data Collection 



In an effort to minimize human errors an established experimental 
procedure was used. The standardized sequence steps were: 

1. Start compressor and wait until operating pressure of 100 psi 
was obtained (15 min.) 

la. Prepare icebath 

2. Adjust air flow rate to achieve desired fluidization 

3. Adjust variac to supply amperage where a temperature 
difference of about 10°C could be achieved 

4. Wait until a steady-state temperature is achieved (10 min.) 
Ha. Check steps 2 and 3 

5. Collect data 

6. Repeat steps 2-5 

Data collection was also performed using a repeated series of events. 
To facilitate data collection and subsequent data reduction, a form 
depicting needed information was used. The form is presented in 
Figure 7. The form indicates locations to record voltage readings for 
the cylinder, bed, icebath and heater cartridge. The dial number is 
the corresponding terminal station on the radial switchbox. It is 
important to note the physical separation of the terminal stations. 
Separation prevents confusion and eliminates the possibility of 
erroneously recording data. The heater cartridge requiring AC volts 
to be read is separated from the thermocouples which require DC volts. 
The thermocouples are separated further into groups pertaining to what 
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temperature they are indicating. A reading is taken by selecting the 
dial number with appropriate voltage mode and observing the voltage 
displayed. This value is then recorded on the form. Fluctuations 
frequently occur because of the violent actions of the fluidized beads. 
When fluctuations occur, the changing values are observed and mentally 
averaged in an effort to isolate the true value. Errors of mental 
averaging are minimized by repeating the readings for five trials. 

Data Reduction 



For a given particle size at each flow rate, a separate data form 
was prepared. Examining Figure 7 reveals locations to reduce the data 
into meaningful information governing the heated cylinder's behavior 
in the fluid. Reduction begins in the upper section by numerically 
averaging the readings at a given dial number for all five trials. 

The subsequent results are then averaged for the cylinder and bed, 
respectively. Using the icebath as a reference, the final bed and 
cylinder voltage values are calculated. This calculated value is then 
converted to a temperature by using the NBS Type K conversion table 
and linear interpolation. The temperature gradient is found by 
finding the difference between the cylinder and bed temperature. 

The heat supplied to produce the temperature gradient is found 
using the electrical expression. 



Q 



v£ 

p 

Q = heat supplied 



where: V = measured cartridge voltage 
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R = known resistance (53-^5 ohms) 

Using the determined temperature gradient, heat supplied and knowing 
the surface area of the cylinder; the desired heat transfer 
coefficient may be calculated. The governing relation is: 



h = 


Q/A 






AT 




where : h = 


overall heat transfer 


coefficient 


A * 


cylinder surface area 


(0.0142 m 2 ) 


AT = 


temperature gradient 


(as measured) 


Q = 


heat supplied 





Nondimensional information is obtained by further data reduction to 
produce the Nusselt number based on particle diameter. 

The Nusselt number based on particle diameter is determined by: 



NU p # = 



.,vp -r 


K 




where: NUp # 


Nusselt number based 


on particle 




diameter 




h = 


overall heat transfer 


■ coefficient 


d P = 


particle diameter in 


m 


K = 


thermal conductivity 


of copper 




(2.63 x 10~ 2 W/m* °C) 





The lower section of Figure 7 lends itself to programmed calculation. 
The large volume of data virtually necessitates programming. 
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To eliminate scale differences of the Nusselt number produced by 
using the various particle diameters as the characteristic dimension, 
the Nusselt number based on the outer diameter of the tube was 
determined. In this manner, a true indication of the total heat 
transfer coefficient is revealed for each particle size examined. The 
constant dimension allows direct comparison of Nusselt numbers between 
data sets. The Nusselt number based on the tube wall diameter was 
calculated by multiplying the existing Nusselt number based on 
particle diameter by the ratio of the wall tube diameter to the 
particle diameter. The governing relation is: 



NU wt - NU p (^) 
a P 



where: NU w t 



NU p 



D'p 



d 



P 



Nusselt number based on tube 
diameter 

Nusselt number based on particle 
diameter for each reading 
outer wall diameter of the tube 
(0.0254 m) 

individual particle diameter in 
meters 



Air flow rates were standardized by determining the mass flow 
rate per minute through the fluidized area. This flow rate, called 
"G" is in units kg/m^S. Flow readings from the three meters were 
converted to the appropriate units using procedures particular to the 



meter used to measure the flow. 



26 



All conversion processes involve knowing the crossectional area 
of the containment vessel. The area is determined by using the well 
known relationship for a circle. Based on a diameter of 11 inches and 
converting the result to m^, the area is 0.063 m^. This value becomes 
a constant for all calculations to determine "G". 

The flow rate for all readings taken on the Small Meter is 
determined with the following relationship: 



G = 

A 



where: % 

C 

P 

A 

M 



= the decimal equivalent of flow read 
directly off the meter 
= capacity of the meter (27 ft3/min.) 

= air density at atmospheric pressure 
and room temperature (1.201 kg/m3) 

» crossectional area of the containment 
vessel (0.025*1 m) 

= metric conversion factor of ft3/min 
to m3/s (0.5886) 



The process governing standardizing readings from the Large Meter is a 
two-step method. The first step involves using a prepared conversion 
curve which permits the rapid conversion of the meter's units into 
SCFM. Briefly, the curve is a plot of percent of flow read off the 
meter versus SCFM. The conversion curve is presented in Figure 8. 

The figure is used by entering with the % of flow read from the meter 
during experimentation on the x axis and recording the corresponding 
flow velocity in SCFM off the y axis. This value is then used in the 
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second step which is a direct calculation. The calculation is 
performed using the relation: 



G 



KpM 

A 



where: K 



P 



A 



M 



value determined from Figure 8 in 
SCFM 

air density at atmospheric pressure 
and room temperature (1.201 kg/m3) 
crossectional area of the containment 
vessel (0.0254 m) 

metric conversion factor of ft3/min 
to m3/S (0.5886) 



Flow rate determination for the Blue Meter is a simple calculation 
involving: 



G 



S&CpM 

A 



where: % 

C = 
P “ 

A = 

M = 



decimal equivalent of flow read 
directly off the meter 
capacity of the meter (1005 ft3/min) 
air density at atmospheric pressure 
and room temperature (1.201 kg/m3) 
crossectional area of the containment 
vessel (0.0254 m) 

metric conversion factor of ft3/min 



to m3/s (0.5886) 
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Using the above procedures results in flow measurements having the 
uniform dimensions of kg/m 2 *S. 

The ratio G/G ra i n forms the dimensionless quantity which is used 
to indicate the relative magnitude of the flow rates used. The G m i n 
value is the smallest flow rate observed to fluidize the specific 
particle being studied. It then becomes the reference point from 
which to measure all other flow rates for that particle. Therefore, a 
G/G m in ratio of 3 indicates the flow rate being used is three times 
the flow rate at the onset of fluidization. 



FIGURE 7 : Standardized Data Collection and Reduction Form 
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RESULTS AND DISCUSSION 



Results 



Data was collected for each particle size at numerous flow rates. 
The assembled results relating the: air flow rate (G Ibm/min* f t 2 ) , 

the ratio of a chosen flow rate to the minimum fluidizing flow rate 
(G/G m i n )> the corresponding overall heat transfer coefficient (h 
W/m 2 '°C) and the dimensionless Nusselt number based on particle 
diameter (NUp) and Nusselt number based on wall tube diameter (NU W ^) 
are presented by particle size. Tables 2-6 present the information 
for particle sizes P-010, P-0230, P-0470, A-205, A-340, respectively. 
Individual result summaries allow trends to be detected for each 
particle size. Relationships which link the results together are 
readily noticed from graphical representations. Figures 9, 10 and 11 
plot results for all particle sizes together. Figure 9 presents the 
relationship of h to G/G m i n . Figure 10 presents Nu p as a function of 
G/G m in* Figure 11 shows the relationship of h to G. Figures 1 3~ 1 6 
plot NU wt vs. G by individual particle sizes. Also indicated on each 
graph, are the predicted results obtained from published correlations. 
The Nusselt number based on tube wall diameter is graphically compiled 



on Figure 12. 
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TABLE 2: Results Summarized for Particle Size P-010 



G min = 0.045 kg/m 2 -S 
d = 2.15 x 10 -4 m 
k = 2.63 x 10"3 W/m°C 



G 


G / G min 


h 


NU p 


NUwt 


0.045 


1 .000 


130.032 


1 .063 


125.582 


0.062 


1.388 


193-741 


1.584 


187.133 


0.075 


1 „ 665 


248.656 


2.033 


240.178 


0.085 


1 .888 


250.647 


2.049 


242.068 


0.101 


2.248 


292.851 


2.394 


282.826 


0.110 


2.442 


306.088 


2.502 


295.585 


0.120 


2 . 665 


309.111 


2.527 


298.539 


0.137 


3.053 


326.214 


2.667 


315.078 


0.160 


3.553 


347.681 


2.842 


335.753 


0.170 


3.775 


351.797 


2.876 


339.769 


0.180 


3.996 


364.600 


2.981 


352.174 


0.194 


4.330 


377.377 


3-085 


364.460 


0.199 


4.440 


383.948 


3.139 


370.840 


0.209 


4.663 


397.218 


3.247 


383.599 


0.219 


4.884 


400.809 


3.277 


387.143 



33 



TABLE 3: Results Summary for Particle Size 



G m in = 0.256 kg/m 2, S 
d = 5.125 x 10 -4 m 
k = 2.63 x 10~ 2 W/m°C 



G 


G/Gmin 


0.256 


1 .000 


0.315 


1 .228 


0.351 


1.370 


0.402 


1 .570 


0.449 


*1 . 750 


0.469 


1 .829 


0.474 


1 .848 


0.517 


1 .992 


0.530 


2.068 


0.564 


2.202 


0.609 


2.378 


0.693 


2.702 


0.785 


3.063 



h 


NU p 


79 . 680 


1 .553 


168.347 


3.281 


219.596 


4.279 


266.210 


5.188 


272.246 


5.305 


273.322 


5.326 


277.788 


5.413 


280.041 


5.457 


279.373 


5.444 


282.021 


5.496 


275.595 


5.370 


273.368 


5.327 


275.730 


5.373 



P-0230 



NU wt 

76.968 
162.610 
212.071 
257.122 
262.921 
263.962 
268. 27 1 ! 
270.454 
269.810 
272.387 
266.1 42 
264.01 1 
266.291 
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TABLE 4: Results Summary for Particle Size P-047 



G min = 0*564 kg/m 2 *S 
d =■ 9.45 x 10 _i4 m 
k = 2.63 x 10~ 2 W/m°C 



G 


O/Gmin 


h 


NU p 


NU wt 


0.564 


1 .000 


136.711 


4.912 


132.026 


0.694 


1.204 


217.665 


7.821 


210.215 


0.774 


1.344 


233.424 


8.387 


225.428 


0.827 


1.435 


236.467 


8.497 


228.385 


0.949 


1.647 


221 .256 


7.950 


213-683 


1.009 


1.750 


223.253 


8.022 


215.618 


1.061 


1.841 


216.326 


7.773 


208.925 


1.108 


1.923 


21 4.844 


7.720 


207.501 


1 .126 


1 .953 


220.619 


7.927 


213.064 


1 .154 


2.003 


216.220 


7.769 


208.818 


1.175 


2.039 


218.988 


7.869 


211 .505 


1.201 


2.083 


215.056 


7.727 


207.689 


1 .243 


2.156 


217.666 


7.821 


210.216 


1.302 


2.259 


213-805 


7.682 


206.479 


CO 


2.563 


209.431 


7.525 


202.259 


1 .709 


2.965 


206.457 


7.418 


199.383 
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TABLE 5: Results Summary for Particle Size A-205 



G min - 1-201 kg/m 2 -S 
d = 2.03 x 1 0“3 m 
k = 2.63 x 10" 2 W/m°C 



G 


G/G min 


h 


NU p 


NU wt 


1 .201 


1.000 


151 .826 


11 .719 


146.632 


1 .31 1 


1 .092 


177.623 


13.710 


171.544 


1.478 


1.195 


196.474 


15.165 


189.750 


1.598 


1.330 


192.909 


14.890 


186.308 


1.662 


1.384 


195.224 


15.053 


188.348 


1 .755 


1.461 


187.899 


14.503 


175.836 


1 .847 


1 .538 


185.272 


14.300 


178.926 


1.939 


1.615 


181.972 


14.046 


175.748 


2.032 


1.692 


181 .494 


14.009 


175.285 


2.124 


1.769 


175.123 


13-517 


169.128 


2.216 


1.846 


170.357 


13-149 


164.524 


2.309 


1.922 


170.367 


13-150 


164.537 



1 
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TABLE 6: Results Summary for Particle Size A-340 



G m in = 1 - 662 kg/m 2, S 

d = 3-40 mm = 3.40 x 10~3 m 
k = 2.63 x 10~ 2 W/m°C 



G 


G/Gmin 


h 


NU p 


NU wt 


1.662 


1.000 


166.788 


21 .562 


161 .080 


1 .755 


1.055 


165.664 


21 .417 


159.998 


1 .847 


1.111 


171 .140 


22.125 


165.287 


1.939 


1.167 


181.715 


23.492 


175.499 


2.032 


1.222 


190.338 


24.606 


183.821 


2.124 


1 .278 


191 .537 


24.761 


184.979 


2.216 


1.333 


191.649 


24.776 


185.091 


2.309 


1.389 


191 .337 


24.736 


184.792 


2.355 


1 .417 


194.510 


25.1 46 


187.855 
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Discussion 



Examination of the results supports the achievement of the 
experimental goal, that is, to identify the total heat transfer 
behavior for all particle sizes specific to the bed under study. 
Published values or correlations vary according to the physical 
make-up of the particular bed studied. Physical factors affecting the 
heat transfer properties include: size, type and weave of the 
distributor plate, the size, shape and thermal properties of the 
fluidized particles, moisture content of the air and lastly, the 
dimensions of the bed. The relationship between particle sizes and 
the trends of a given bead size are most readily noticed from 
graphical displays. Figure 11 shows the behavior of the total heat 
transfer coefficient as the absolute flow increases for all bead sizes. 
Figure 10 depicts the characteristic Nusselt number based on particle 
diameter for flow rates above minimum fluidization. Figure 9 shows h 
as a function of relative flow. All particle sizes are plotted 
together. To identify the individual particle sizes, different data 
symbols are used. The symbols are presented in Table 7 and on each 
graph for clarity. 

Examining the plot of h vs. G/G m i n or Figure 9 produces 
interesting observations. The smallest particle sizes tend to have 
the highest total heat transfer coefficient. This behavior is most 
evident at G/G m i n values above 1.4. So, as particle size increases, 
total heat transfer decreases, in agreement with many previous 
investigations [2,4,7]. The decrease of total heat transfer with 
particle size is due to the increasing importance of gas convection 
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for the larger particle sizes. By tracing each particle size 
individually, a repeated pattern is noticed. The pattern tends to be 
one of rapid increase of h for a reasonably small change in G/G m i n 
until a maximum h is achieved at a specific G/G m i n . After this unique 
point, the h value decreases slowly for higher flow rates. Particles 
P-047 and A-205 demonstrate the behavior most dramatically. Bead size 
P-010 was not observed in a flow rate range large enough to include 
its maximum. 

TABLE 7 : Particle Name and Corresponding Data Symbol 



SYMBOL 


BEAD 

IDENTIFIER 


V 


P-010 


o 


P-023 


□ 


P-047 


O 


A-205 


A 


A-340 



The uniqueness of a particular flow rate corresponding to a 
maximum h value is clearly seen in Figure 11. On this figure, the G 
at the largest total heat transfer coefficient may be read directly 
from the graph. The value of fluidizing velocity at which the maximum 
in h occurs is called the optimum velocity, U 0 p£. It has been 
observed by many investigators [10] and a correlation provided by 
Toades [11] may be used to predict its value. The correlation 
presents the optimum Reynolds number. The correlation yields the 
corresponding U 0 pp value to within a multiplied constant involving the 
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specific kinematic fluid viscosity and the characteristic length of 
the apparatus used. The expression becomes: 



» Ar r Y-) 

°P fc 18 + (5) (22/Ar) ''L' 1 

where: Ar = the Archimeder number 

Y * kinematic viscosity of air 
(1.55 x 105 m 2 /S) 

L = characteristic length or the diameter 
of each particle size 

This relationship was applied for all particle sizes and the results 
compared to the experimental results of the present study. Table 8 
compiles the information by particle size. 

TABLE 8 : Summary of U 0 p^ by Particle Size 



U opt (m/s) 



Beads 


Experiment 


Correlation 


P-010 


Not Applicable 


0.3580 


P-023 


0.4225 


0.5972 


P—0‘47 


0.6886 


0.8379 


A-205 


1.2306 


1 .2293 


A-340 


1.6919 


1.5941 



The correlation velocities fit the experimental results very well for 
the largest two particle sizes, however the agreement is not quite as 
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good for the intermediate particle sizes. Total heat transfer 
decreases at velocities above the optimum velocity because of the 
increasing importance of gas convection heat transfer. As the 
fluidizing velocity is increased, a corresponding increase in bed 
porosity occurs. The excess fluid separates the particles which 
diminishes the strong contribution of particle heat conduction, 
thereby lessening the overall heat transfer. 

Closer examination of the h vs. G plot shown in Figure 11 
indicates that the dependency of h on G is stronger for smaller 
particle sizes. That is, the slope of the curve is greater for small 
beads. Using Figure 9, the steep rise is obvious when beads P-023 and 
P-047 are compared to A-205 and A-3^0. For G/G m i n values from 1.0 to 
about 1.6, the h values of the smaller particles changed from about 80 
to 265 and 135 to 236, respectively as compared to 150 to 195 and 1 65 
to 190, respectively for the larger particles. 

This behavior may be explained by considering the mechanisms 
involved which compare the overall heat transfer coefficient. It is 
the sum of a conductive part and gas convective part. For smaller 
particles, the h value is dominated by particle conduction. Heat 
transfer predominantly caused the physical touching of packets of 
particles to the heated cylinder. Whereas for the largest particles, 
the overpowering mechanism is gas convection. 

Figure 10 plots the Nusselt number based on particle diameter 
versus relative flow rate for all particles. The plot clearly 
indicates that each bead size has a characteristic range of Nusselt 
number. The band tends to be narrow and at successively higher values 
as bead size increases. The increase in Nusselt number with 



41 



increasing particle size does not indicate conflicting results as 
presented in Figure 11. The Nusselt number calculation is controlled 
by the particle diameter dimension. This dimension overpowers the h 
value. The general shape of the curves for both figures is consistent 
with the earlier work of others [2,4,7]. 

Nusselt numbers based on particle diameters can be somewhat 
misleading. They may erroneously imply that heat transfer increases 
with increasing bead size. Nusselt numbers based on the tube wall 
diameter accurately portray the physical occurrence since the tube 
wall diameter is constant. Compare the implied behavior of Figure 12 
for NU wt as a function of G using all particle sizes to that of Figure 
10. As expected, the plotted curves have shapes identical to the 
plots of h on G. Figure 12 shows that the total heat transfer 
coefficient increases with decreasing bead size. 

Detailed examination of NU W ^ in relation to G for individual 
particles is presented in Figures 1 3~ 1 7 . Each display also has the 
behavior predicted by the three most applicable correlations. The 
correlations proposed by Grewal and Saxena [7], Andeen and Glicksraan 
[4], and Vreendenberg [2] are shown below. Each investigator 
attempted to account for the many factors which affect heat transfer 



in a fluidized bed. 
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Grewal & Saxena: 



0.325 



. ,7(1 - E )(S^ g^) ' 



p T 3/2 g 1/2 ^ 



23 



k f 



p r 0. 30 



where 



CTto.* 



r 0.43 1/3 

[4[ d?rp7[7s Hr Pf)^ g j } } ] 



Andeen & Glicksman: 



NU wt = 900 (1 



.urGDtgau,, m 2 n 
' 1 PfP 



0.326 

Pr 0.3 



where: 



e 




200 Gy 

Pflps ~ PfJg 



} 



1/3 



Vreendenberg: 

0 3 

NU wt = 420[(^-^)( ^/ - 2 -)] ' Pr0-3 
wc PfU d p? g ; 

The nomenclature common to all the correlations and the values 
which correspond to this experimental setup are as follows: 



e = bulk porosity (as determined above) 

G = superficial mass fluidizing velocity (G in. kg/m^'S) 

D<j* - outside tube diameter (0.0254 m) 

p s = density of the bed material - for glass beads (2500 Kg/m^) 

pf = density of the fluidizing air (1.201 kg/m3) 
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y 


fluidizing air viscosity (184 x 10"? N*s/m 2 


d p = 


particle diameter as shown on TABLE 1 


g 


acceleration due to gravity ( 9 . 8 1 m/s 2 ) 


Cps = 


specific heat of the bed material - for glass beads (750 J/kg-°K) 


Kf = 


thermal conductivity of the fluidizing air (26.3x10”^ W/m-°K) 


Pr 


Prandtl number of air (0.707) 


4>s “ 


sphericity of the beads (assumed to be 1 ) 



For the 


smallest particle size, P-010, the Andeen and Glicksman 



correlation most accurately represents the data. Figure 13 shows the 
experimental NU w p values and the correlation values as a function of G. 
This representation is about 1056 different from measured values. 

Figure 14 presents the same information for P-023* Here, the Grewal 
and Saxena correlation coupled with the Andeen and Glicksman 
correlation form the upper and lower bounds, respectively. Both 
correlations could be interpreted to represent the data. The Grewal 
and Saxena correlation is most accurate for intermediate superficial 
velocities. At larger flow rates, the Andeen and Glicksman model is 
most applicable. It should be noted that none of the correlations 
account for the rapid increase of the heat transfer coefficient at the 
smallest flow rates for either P-010 and P-023. 



Figures 


15, 16, 17 present the experimental and correlated 


behavior for 


bead sizes P-047, A-205 , A-340, respectively. The 


experimental 


results are very poorly predicted by the three 


mathematical 


models in all cases. The correlations anticipate NU wt 



values significantly lower than those actually achieved. The 
correlations were intended to be used for particles less than 500 
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micrometers. Since the last three bead sizes clearly exceed the 
particle size limit, extrapolation may not be valid. The results are 
qualitatively compatible and quantitatively substantiated by the work 
of others in the field. 
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FIGURE 10 s Nusselt Number Based on Particle Diameter Against i 
Relative Mass Velocity 
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FIGURE 11 i Heat Transfer Coefficient Against Absolute Mass Velocity 
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FIGURE 13 : Comparison of Experimental to Emperical 

Nusselt Numbers Against Absolute Mass 
Velocity for P-010 
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FIGURE 14: Comparison of Experimental to Emoirical 

Nusselt Numbers Against Absolute' Mass 
Velocity for P-023 
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FIGURE 17 : Comparison of Experimental to Emperical 

Nusselt Numbers Against Absolute Mass 
Velocity for A-340 
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CONCLUDING REMARKS 

Experiments were conducted to identify the heat transfer 
characteristics particular to the fluidized bed studied. The goal was 
to quantify the total heat transfer coefficient as a function of the 
mass velocity for five different fluidized particles. The particles 
ranged in size from 2.15 x 10 -14 m to 3.40 x 10“3 m . Examination of 
tabulated and graphical data revealed the following generalizations. 

The smallest particle sizes tend to have the highest total heat 
transfer coefficient. So, as bead size increases, the overall heat 
transfer decreases. The decrease of total heat transfer with 
increasing particle size is attributed to the successively growing 
importance of gas convection. 

Each particle size exhibited a similar h versus G pattern. The 
pattern is characterized by a rapid increase in h for a small change 
in G until a maximum h is achieved at a specific G; followed by a 
slowly diminishing h value over all remaining G values. The patterns 
first stage clearly demonstrates a stronger dependency of h on G for 
smaller fluidized bead sizes. The superficial velocity which 
corresponds to the maximum heat transfer coefficient in each case is 
the optimum velocity for that particle size. This optimum velocity, 
U 0 pt waa experimentally identified for each bead size. Knowing this 
ideal point will aid subsequent heat transfer experimentation using 
the same apparatus. The experimentally determined U 0 pt values compare 
favorably to empirically defined values found from the Todes 
correlation [11]. The last section of the pattern shows h slowly 
tailing off. The gradual decrease is caused by the increasing bed 
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porosity at higher velocities. The larger fluid quantities decrease 
particle heat conduction, thereby enabling gas convection to dominate 
the heat transfer process. 

The dimensionless Nusselt number was used as a basis for analysis. 
The Nusselt number based on particle diameter, NUp/f, shows that each 
bead size has a narrow range of NU p # values with the values increasing 
with particle size. The Nusselt number based on tube wall diameter, 
NU w t#> was used to directly compare the heat transfer data at all 
velocities among the various particle sizes. NU W ^# behaved in a 
manner identical to the behavior of h. The validity of the NU w t# 
values was tested against results obtained from numerical correlations 
by Grewal and Saxena, Andeen and Glicksman, and Vreendenberg. The 
correlations substantiated findings for the smallest three particles. 
For the largest two particles, the correlations do not reflect the 
effects of the large beads. This makes them invalid and were not 



applied in the analysis. 
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